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We numerically investigate Andreev reflection in a graphene ring with one normal conducting 
and one superconducting lead by solving the Bogoliubov-de Gennes equation within the Landauer- 
Biittiker formalism. By tuning chemical potential and bias voltage, it is possible to switch between 
regimes where electron and hole originate from the same band (retroconfiguration) or from diflerent 
bands (specular configuration) of the graphene dispersion, respectively. We find that the dominant 
contributions to the Aharonov-Bohm conductance oscillations in the subgap transport are of period 
hj2e in retroconfiguration and of period h/e in specular configuration, confirming the predictions 
obtained from a qualitative analysis of interfering scattering paths. Because of the robustness against 
disorder and moderate changes to the system, this provides a clear signature to distinguish both 
types of Andreev reflection processes in graphene. 

PACS numbers: 72.80.Vp, 73.23.-b, 74.45. +c, 85.35.Ds 



Since its first experimental realization in 2004 [T], 
graphene has strongly influenced the field of mesoscopic 
physics due to the peculiar behavior of its electronic ex- 
citations as massless Dirac fermions. Two prominent 
features, namely the effects of Klein tunneling [2] and 
Andreev reflection [3], concern processes at interfaces 
between regions of different doping and normal metal- 
superconductor (NS) junctions, respectively (for a re- 
view see Ref. While both effects are known to 
be closely related [S], different aspects of Klein tunnel- 
ing have become experimentally accessible in the last 
years [^, whereas specular Andreev reflection has not 
been observed to date, although there exist a number 
of proposals for the experimental control [7] and detec- 
tion [21 [S] of this process. In this Letter, we will present a 
novel approach concerning the identification of specular 
Andreev reflection, distinguishing it from conventional 
retroreflection, and discuss the advantages over previous 
works in the field. 

Our approach is based on the observation that in gen- 
eral, the probability for an incident electron to be re- 
flected as a hole is less than one. This allows for effects 
typical for phase-coherent mesoscopic devices, like uni- 
versal conductance fluctuations or Aharonov-Bohm os- 
cillations in the magnetoconductance. While in nor- 
mal metals, the fundamental period of these oscillations 
is given by the flux quantum $o = h/e, it is half the 
value for Andreev (retro)reflection in conventional met- 
als, due to the charge 2e of a Cooper pair. However, 
this is not true anymore in the case of specular Andreev 
reflection, therefore providing a criterion to distinguish 
between specular and retroreflection. In order to show 
this, we consider the phases due to the magnetic flux 
that are picked up by the various scattering paths. In 
this analysis, we restrict ourselves to the contributions 



up to first order in the sense that we take processes into 
account that involve only a single electron-hole conver- 
sion process, and that contain at most one additional 
round-trip of electron or hole, respectively; higher order 
contributions connected with additional round-trips are 
often times negligible [lOl [11] . The corresponding paths 
are summarized in Fig. [l] In order to obtain the mag- 
netoconductance for the two types of Andreev reflection 
[specular (s) and retro (r)], we sum up the amplitudes as 
defined in Fig. [l] for the various paths coherently to ob- 
tain the corresponding Andreev reflection probabilities: 

i?,($) - + .s_ + s'+c'* + s'_e-'*f 

i?,($) ^ |r+e'* + r_e-'* + r^e^'* + rLe-2'*|' (1) 

where s'j_ — s'j_g + s'^f^, = r^g -I- r^^, and is the 
magnetic flux measured in units of the flux quantum $o- 
Assuming \s\ 3> \s'\ for any zeroth- and first-order am- 
plitudes, respectively, we obtain 

i?,($) - i?o + 2Re [{s'^sl + sos'!)e'*] + 0[{s')\ (2) 

where sq = ,s+ + .s_ and i?2 contains contributions that 
are constant with respect to $. Therefore, in the case 
of specular reflection, oscillations of period h/e are dom- 
inant. In contrast, in the case of retroreflection, contri- 
butions of period h/2e are dominant, as expected: 

i?,($) ^ i?° + 2Rc [r+rlc^'*] + 0[rr' , {r')\ (3) 

where again i?" contains ^-independent terms and we as- 
sume |r| ^ \r'\ for any zeroth- and first-order amplitudes, 
respectively. 

In order to test this analysis on the basis of a micro- 
scopic model, we implement the Bogoliubov-de Gennes 
Hamiltonian [T^ 
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FIG. 1. (a) Device geometry showing a graphene ring struc- 
ture that is penetrated by a magnetic flux <E> measured in 
units of the flux quantum $o- At the interface with the su- 
perconductor (shaded region), electron- hole conversion may 
occur, (b) The gauge is chosen such that each of the eight 
individual electron (solid lines) and hole (dashed lines) paths 
picks up a phase ±cl?/2 as indicated, (c) Scattering paths for 
electrons injected from and holes leaving through the left nor- 
mal conducting lead; only zeroth and first order contributions 
are included, i. e. terms containing a single electron-hole con- 
version process and at most one additional round-trip of the 
electron or the hole. The paths are categorized according to 
the total phase that is picked up, and each path is associated 
with a corresponding amplitude, where first order amplitudes 
are indicated by a prime. 



within the tight binding formalism of graphene 



(5) 



where the second sum runs over nearest neighbors and 
Ui = U{Yi) is a position-dependent potential. In this 
numerical calculation, all higher order contributions are 
also taken into account. In Eq. ([4|, we assume A, = 
A(ri) € M for the superconducting order parameter A = 
^ ■ Ai |i) The presence of a magnetic field is captured 
by a Peierl's phase in the hopping matrix element 



-To exp 
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A(r)dr 



(6) 



where Tq ~ 2.7 eV is the graphene hopping integral, $o = 
h/e is the magnetic flux quantum, and the line integral 
is taken along the straight path between sites i and j. 

The structure of the graphene device under consider- 
ation is schematically shown in Fig. [l] The ring-shaped 
structure is generated by setting the appropriate hop- 
ping matrix elements to zero in Eq. The two semi- 
infinite leads also exhibit the graphene lattice structure; 



superconductivity is induced into the right lead due to 
the proximity effect of a superconducting electrode on 
top of the graphene. We choose to orient the leads to ex- 
hibit armchair edges and later comment on the reason for 
this particular choice. The whole ring is penetrated by a 
uniform perpendicular magnetic field of strength _B, de- 
scribed by the vector potential A(r) = —By9{d — \x\)ex- 
The origin of coordinates is taken at the center of the 
ring. 

In order to fulfill the mean-field requirement of super- 
conductivity, which demands the superconducting coher- 
ence length ^ = hvp /A to be large compared to the wave- 
length \s in the superconducting region [3] , we introduce 
additional doping into the superconducting region by ap- 
plying a gate potential Ui = U9{xi — d). Which type of 
Andreev reflection occurs at the NS interface is then de- 
termined by the excitation energy e [i. e. the eigenvalues 
of Eq. (|4])] and the Fermi energy Ep, as shown in Fig. [2] 



In retroconfiguration, Ep > e > 0, where Vy 
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Uy 
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for the y-components of the electron and hole velocities, 
both electron and hole traverse the same arm of the ring. 
In specular configuration, < Ep < e, the hole is re- 
flected back through the other arm of the ring, since 
Vy'^'^ ■ v'^^ > 0. In the following, we choose \U\ ^ Ep, 
justifying the adoption of the step- function model for the 
superconducting order parameter, A^ = A0(xi — d) |5- 

In order to compare retro (r) and specular (s) config- 
urations, we will choose e^^^ = E^'p' and e'"-' = i?^'* since 
then the states in both configurations exhibit the same 
wavelength and there is the same number of propagat- 



ing modes. We further choose £^'^\Ep ^ £^^\Ep so 
that for nearly each value of ky , there exist electron- hole 
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FIG. 2. Schematics of the excitation spectrum (lower panel) 
and surfaces of constant excitation energy in k space (up- 
per panel) in the cases Ep > e > (retroconfiguration) and 
< Ep < e (specular configuration). Solid and dashed lines 
indicate electron- and hole-like states, respectively, (hole) 
states originating from the valence band are shaded gray. The 
small arrows in the upper panel indicate the direction of prop- 
agation of the corresponding states. Electron-hole excitations 
are drawn assuming conservation of ky at the NS interface. 
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scattering channels. 

The transport properties of the system are obtained 
from the scattering matrix S that is calculated in the 
framework of the Landauer-Biittiker formalism using a 
variant of the recursive Green's function technique (13] 
we recently applied to a similar setup [10] . This technique 
is an efficient way to obtain the relevant parts of the sys- 
tem's Green's function from the surface Green's functions 
of the isolated leads, which is known analytically [H], by 
solving Dyson's equation exactly, treating the coupling of 
the leads to the ring region as perturbation. The Fisher- 
Lee relation [15] then relates the Green's function to the 
S matrix, from which the transmission function may be 
obtained. In this framework for elastic transport. Green's 
function and scattering matrix are parameterized by the 
eigenvalues e of the Hamiltonian Q. 

In the following, we will concentrate on the regime 
e < A, in which there are no propagating modes in the 
superconducting lead, so that electrons injected from the 
normal conducting lead are reflected back either as elec- 
tron (e) or hole {h). The scattering matrix thus has the 
structure 



The rhh 



(7) 



from which the differential conductance for the Andreev 
processes is given by 



dl 
dV 



4e2 



(8) 



where the factor 4 accounts for spin degeneracy and the 
quantization of charge in units of 2e. 

In Fig. |3] we show the calculated transmission for a 
ring of width w = 87-\/3ao and outer radius R = 500ao, 
where ag is the distance between nearest neighbors. The 
transmission function exhibits Aharonov-Bohm oscilla- 
tions on top of a low frequency background which is due 
to universal conductance fluctuations. The position of 
the NS interface is given by d = 400ao. The chosen 
dimensions of the ring are large enough to exclude finite- 
size effects while still being numerically manageable. For 
the superconducting order parameter, we choose a value 
of A = O.OStq « 80 meV, which may appear unrealistic 
at first sight, considered the fact that typical values are 
up to a few meV. However, by making this choice we scale 
the value of the superconducting order parameter accord- 
ing to the scale of the system size, such that the dimen- 
sionless factor AR/Hvp stays of same order of magnitude, 
compared with values realized in experiments |111 I16j . 
Thus, for a realistic system size of i? ^ lO^^m, our choice 
of A would correspond to a value of a few meV for the 
superconducting gap. Note that due to these low energy 
scales and the rather large spacing of modes resulting 
from the narrow geometry of the electron waveguides in 
such a ring structure, in specular configuration only the 
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FIG. 3. Differential magnetoconductance for specular (black) 



and retro (gray) configuration for Ep 
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Ep — O.OOIto — e^''\ corresponding to 8 modes in the nor- 
mal conducting lead, including all degeneracies (spin, valley, 
electron or hole). The high doping in the superconducting 
lead is chosen such that Ef — U = 0.5to in both cases. Other 
parameter values are provided in the main text. The period of 



the dominant oscillation is B, 
ular configuration and B^'' 
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in retroconfiguration. The weak beating pattern in retrocon- 
figuration and the asymmetry in specular configuration arise 
due to minor contributions of contrary frequencies. 



regime of a low number of modes is accessible for realis- 
tic choices of system parameters. Also note that due to 
strong electron backscattering at the front of the hole and 
at the rough edges of the ring, the average value of the 
differential conductance is much less than a conductance 
quantum, e'^/h. 

The average radius f of the scattering path is calcu- 
lated according to t^tiBq = h/ne, where n = 1 (n = 2) 
in specular (retro) configuration and Bq is the (domi- 
nant) period of the oscillation. Evaluating the period of 
the oscillations shown in Fig. la) we obtain r^^) « 420ao 
in specular configuration and f ~ 425ao in retrocon- 
figuration. The obtained values lie well within the inner 
and outer radius of the ring and close to the arithmetic 
mean R — w/2 « 425ao, therefore confirming the pre- 
dictions obtained from Eqs. ([2| and (|3]). Minor contri- 
butions of period h/e in retroconfiguration and h/2e in 
specular configuration visible in Fig. [3] may arise due to 
terms neglected in Eqs. ^ and ([s]), scattering off the 
sharp boundaries of the ring structure, and the fact that 
for the electron-hole conversion at the NS interface ky is 
not strictly conserved. 

Other strong evidence that supports our interpreta- 
tion of the two different periods is the breakdown of this 
particular signature that is observed for a shift of the po- 
sition of the NS interface on the scale of the width of the 
ring. Indeed, while in Ref. [17] — where a three-terminal 
graphene junction is analyzed — the exact position of the 
NS interface has no effect, it matters in our case; the 
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FIG. 4. Breakdown of the h/e vs. h/2e signature for shifted 
positions of the NS interface, as explained in the text. Other 
parameters and color coding are chosen as in Fig. [3] For 
d = 340ao (left), in specular configuration one observes oscil- 
lations of period h/2e as in retroconfiguration. For d = 490ao 
(right), contributions of specularly reflected holes in retro- 
configuration become important, leading to the observation 
of additional /i/e-oscillations. The value of the superconduct- 
ing coherence length is ^ = SOao- 



reason is that ^ is comparable or even less than the sys- 
tem size, while in Ref. [17^ the superconducting coherence 
length greatly exceeds the system dimensions. If the in- 
terface is too close to the hole region [see Fig. |4]^a) inset], 
then specularly reflected holes are forced to traverse the 
same arm as the incoming electron. In this case, one 
should observe h/2e oscillations in specular configura- 
tion. In contrast, if the interface is too far from the hole 
[see Fig. |4|^b) inset] , holes may significantly be reflected 
through the other arm, e. g. due to increased scattering 
at the ring boundaries. This would manifest itself in the 
observation of h/e oscillations in addition to the h/2e 
oscillations in retroconfiguration. This behavior is con- 
firmed in the observed magnetooscillations, as shown in 
Fig. [4] 

Apart from that, the h/e vs. h/2e signature proves to 
be very robust against moderate changes to the length 
and energy scales in the system, such as the extent of the 
magnetic field or the ratio of Fermi wavelength and the 
width of the NS interface. We also tested that the signa- 
ture persists when more propagating modes are present 
in the lead, leading to values of the average conductance 
which are much larger as compared to the few-mode sit- 
uation shown in Fig. |3] Additionally, the signature is 
hardly affected by bulk disorder, which is a major ad- 
vantage of our setup. In Fig. [5| we show the magne- 
toconductance of the system used in Fig. [3] with a par- 
ticular random short-range disorder configuration, which 
is realized by applying an uncorrelated, random on-site 
potential of Gaussian distribution with zero mean and 
width a = O.OItq to each site. In addition, the NS in- 
terface has been smeared out over a distance I — 90ao in 
this case. The robustness of the effect can be explained 
from the topological nature of the signature: since all mi- 
croscopic scattering paths can be classified into just two 
groups — yielding h/e- or ft,/2e-oscillations, respectively — 
according to which arm is traversed by the quasiparti- 
cles, impurity scattering and the resulting defiection of 



quasiparticles has no adverse effect as long as scatter- 
ing between the groups is weak, while scattering within 
one group may be arbitrarily strong. In addition, note 
that while our description of transport via the scattering 
matrix assumes complete phase coherence, a signature 
that distinguishes retro from specular Andreev reflection 
is assumed to persist also in the case of a flnite phase co- 
herence length. More specifically, if the phase coherence 
length is on the order of the ring circumference, first- 
order amplitudes in Eqs. ^ and (|3| may be neglected. 
Then, retrorefiection would still manifest itself in h/2e- 
oscillations, while there would be no oscillations at all in 
the case of specular reflection. 

Before we conclude, we like to add a remark concern- 
ing the choice of armchair boundary conditions in the 
leads we employed in our analysis. In a tight binding im- 
plementation of graphene, there are two simple choices 
for the orientation of the leads. Often, zigzag edges 
are considered to represent a generic boundary condi- 
tion for graphene ribbons [TS]. In this case, edge states 
are present in the system that modify the simple pic- 
ture provided in Fig. [2] by adding additional scattering 
channels between bulk and edge states while removing 
certain scattering channels between bulk states due to 
the conservation of the so-called pseudoparity symmetry 
that acts like a selection rule [H] . In the realistic limit of 
metal leads providing a large number of propagating bulk 
modes, this effect should be less important. However, for 
the system geometry used in the numerical calculations 
in combination with the low energy scales, it may signif- 
icantly affect the observed behavior. In order to avoid 
this influence, we chose armchair boundary conditions in 
the leads that do not provide any edge states. Note in ad- 
dition, that in realistic systems the zigzag-speciflc effect 
would also be suppressed since the zigzag edge state is 
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FIG. 5. Magnetoconductance of the system used in Fig. [3] 
with a smooth potential profile (inset) with / — 90ao and 
bulk disorder of strength a = O.OItq as explained in the text. 
The /i/e vs. h/2e signature still persists. The color coding is 
the same as in Fig. |3] 
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not protected against disorder when next-nearest neigh- 
bor hopping is taken into account [20j . Therefore, we are 
convinced that our results based on armchair edges in the 
reservoirs describe the generic situation for wide leads. 

In conclusion, we have shown numerically that the fre- 
quency of Aharonov-Bohm oscillations in graphene rings 
provides a clear and feasible signature for distinguish- 
ing specular Andreev reflection from retroreflection. This 
feature can be explained qualitatively by considering the 
interference of the different scattering processes up to 
first order. The signature is robust against the presence 
of disorder and persists within a certain range for the 
position of the NS interface before it breaks down when 
the interface gets too close to or too far away from the 
hole of the ring. 
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